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ABSTRACT 

We have carried out observations of the X-ray transient GX 339—4 during its high-soft 
and low-hard X-ray spectral states. Our high-resolution spectroscopic observation in 
1999 April suggests that the Ha line has a single-peaked profile in the low-hard state as 
speculated in our previous paper. The He ii A 4686 line, however, has a double-peaked 
profile in both the high-soft and low- hard states. This suggests that the line-emission 
' mechanism is different in the two states. Our interpretation is that double-peaked lines 

\ are emitted from a temperature-inversion layer on the accretion-disk surface when 

. it is irradiatively heated by soft X-rays. Single-peaked lines may be emitted from 

outflow/wind matter driven by hard X-ray heating. We have constructed a simple 
plane-parallel model and we use it to illustrate that a temperature-inversion layer can 
be formed at the disk surface under X-ray illumination. We also discuss the conditions 
O ' required for the formation of temperature inversion and line emission. Based on the 

velocity separations measured for the double-peaked lines in the high-soft state, we 
^2 ' propose that GX 339—4 is a low-inclination binary system. The orbital inclination is 
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about 15° if the orbital period is 14.8 hours. 

Key words: binaries: spectroscopic — stars individuals (GX 339—4) — black hole 
physics — accretion: accretion disks 



1 INTRODUCTION 

GX 339—4 is an X-ray transient (Markert etal. 1973; Harmon etal. 1994), and is classified as a black hole candidate (BHC) 
because of its short-term X-ray and optical variability (Makishima et al. 1986) , its transitions between high-soft and low-hard 
X-ray spectral states (Markert et al. 1973) , and the extended high-energy power-law tail in its X-ray spectrum (Rubin et al. 
1998). It is one of the few BHCs that has shown all four X-ray spectral states: ojj, low-hard, high-soft and ultra-high. During 
the off state, GX 339—4 has very weak and hard X-ray emission. Its optical counterpart is faint, with 1/ ~ 19 — 21. In the 
low-hard state, it emits very strong hard X-rays. The spectrum is an extended power-law with a photon index ~ 1.5. The 
2—10 keV X-ray flux is, however, weak (~ 0.4 x 10^^ ergcm"^ s~^). The optical brightness of the system is 1/ ~ 16 — 17. In 
the high-soft state, the X-ray spectrum is characterised by a soft, thermal component. The 2—10 keV X-ray flux is about 20 
times higher than in the low-hard state. The power-law tail is weak, with a photon index ~ 2. The optical brightness is similar 
to that in the low-hard state, with ~ 16 — 17. In the ultra-high state, both the thermal and the power-law components in 
the X-ray spectrum are strong. The 2—10 keV X-ray flux is about 50 times the X-ray flux in the low state, and the photon 
index of the power law is ~ 2.5. On one occasion Mendez & van der Klis (1997) reported that the source was in a state 
intermediate between the low-hard and the high-soft state. The 2—10 keV X-ray flux was about a factor of 5 below the X-ray 
flux in the high state. 

GX 339—4 is also a radio source (e.g. Hannikainen etal. 1998). It has a flat spectrum, and its image has shown jet-like 
features (Fender etal. 1997). The radio emission was found to be correlated with the hard (20 — 100 keV) X-ray flux, and was 
quenched during a high-soft state (Fender et al. 1999) . 

Despite the fact that GX 339—4 is optically bright [V ^ 1& ~ 17) when it is X-ray active, it is not well studied in the 
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Table 1. Log of our oj)tical sj)Cftroscoyjic oljsorvatioiis. 



UT date 


epoch of observations 


Wavelength range 


Resolution 


Telescope 




(HJD - 2450000) 


(A) 


(A FWHM) 




low-hard X-ray spectral state 


1997 May 6 


574.966-575.226 


5355-6950 


3 


AAT 


1997 May 8 


576.974-577.042 


5355-6950 


3 


AAT 


1999 April 12 


1281.287-1281.313 


4455-5072 


1.3 


AAT 






6282-6828 


1.3 






liigii-sui'l X-ra\' sjx' 


(iral si ale 






1998 April 28 


931.953-932.318 


4150-5115 


1.3 


ANU 2.3 m 






6200-7150 


1.3 




1998 April 29 


932.943-933.268 


4150-5115 


1.3 


ANU 2.3 m 






6200-7150 


1.3 




1998 April 30 


934.106-934.233 


4150-5115 


1.3 


ANU 2.3 m 






6200-7150 


1.3 




1998 August 20 


1045.864-1046.140 


4150-5115 


1.3 


ANU 2.3 m 






6200-7150 


1.3 




1998 August 23 


1048.860-1049.103 


4150-5115 


2 


ANU 2.3 m 






6200-7150 


2 





optical bands. The mass function of the system has not yet been determined, and so the black hole candidacy is not verified in 
terms of the orbital dynamics. In a photometric observation during the off state, Callanan etal. (1992) detected a periodicity 
of 14.8 hr, and attributed it to the orbital period. A 14.8-hr orbital period would imply that the system has a late-type 
companion star with a mass < 1.6 Mq. The 14.8-hr period, however, was not seen in the photometric data that we collected 
in 1998 (Soria 2000; Soria, Wu & Johnston 2000). 

In our previous paper (Soria, Wu & Johnston 1999, hereafter Paper I) we reported part of the findings from our 1997 
and 1998 spectroscopic observations. The system was in a low-hard state during the 1997 observation and was in a high-soft 
state during the 1998 observations. We found that the emission lines generally show double-peaked profiles in the high-soft 
state. The Ha line did not show resolved peaks in the medium-resolution (3 A) spectra that we obtained in the low-hard 
state. Because of the limited spectral coverage in our 1997 observation, we were unable to determine the properties of the H/3, 
Hen A 4686 and Bowen Niii AA 4641,4642 lines in the low-hard state. We therefore restricted our discussion to the profiles of 
the Ha line and did not compare the properties of the high-ionisation lines in the two X-ray spectral states. 

GX 339—4 was in its low-hard state again in early 1999, but it entered an X-ray quiescent phase in mid- April, with its 
hard X-ray flux below the detection limit of BATSE (from BATSE bright source reports) and its RXTE/ASM count rate 
consistent with zero (from quick-look results of the RXTE/ASM team). Its radio flux density in late June, measured by the 
ATCA, was more than 10 times weaker than that measured in an observation on April 14 (S. Corbel, private communication). 
Our optical spectroscopic observation was carried out on April 12, about 5 days before it became undetectable by BATSE. 
High- resolution (~ 1 A) spectra at the Ha and the H/3/Heii regions were obtained. In this paper we report the results of 
the 1999 observation and compare them with those of our 1997 and 1998 observations. We discuss the profiles of the high- 
and low-ionisation lines in different X-ray spectral states and the implications for the line-formation mechanisms and the 
line-emission region in BHCs. 



2 OPTICAL SPECTROSCOPY 

Our spectroscopic observations were carried out with the 3.9-m Anglo-Australian Telescope (AAT) in 1997 May and 1999 
April and with the ANU 2.3-m Telescope at Siding Spring Observatory in 1998 April and August. The observations covered 
two X-ray spectral states of the source: the low-hard state in 1997 and 1999, and the high-soft state in 1998. 

The observations are listed in Table 1. The technical details of the 1997 and 1998 observations were presented in Paper I. 
For the 1999 April observations, we used a Tektronix Ikxlk CCD on the RGO spectrograph and a 1200V grating; the seeing 
was < 1 arcsec. The spectra were centred at 6563 A (Ha region) and at 4770 A (H/3/Heii region), and the spectral resolution 
was 1.3 A FWHM for both the Ha and the H/3/Heii spectra. The total exposure time for the Ha spectrum was 1200 s, and 
that for the H/3/Heii spectrum 600 s. Standard data reduction procedures were followed, using IRAF tasks, and wavelengths 
were calibrated using the Cu-Ar lamp spectra. 
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Figure 1. The profiles of the Ha, H/3, Hei A6678, Hei A7065, Hen A4686 and Niii AA 4641,4642 fines of GX 339-4 in 1998 April 
during a high-soft state. 

3 LINE PEAK-SEPARATION IN THE HIGH-SOFT STATE 

In Paper I wc showed the Ha and H/3/Hen spectra of GX 339—4 in 1998 April and August. The strong emission lines in the 
spectra are Hi Balmer, Hei A 6678, Hei A 7065, Hen A 4686 and Bowen Niii AA 4641,4642. Other detectable emission lines are 
Hei A 4922, Hei A 5876 and Nii A 6505. The spectra are very similar to those of other BHCs (e.g. GRO J1655-40, see Soria, 
Wu & Hunstead 2000) in their high-soft states. The spectra of these systems all show high-ionisation lines, and the emission 
lines are double-peaked. 

We show in Figure 1 the high-resolution profiles of the double-peaked Ha, H/3, He: A 6678, Hci A 7065, Hen A 4686 and 
Niii AA 4641,4642 lines. The velocity separations of the peaks are 250±27 km s"^ for Ha, 260±35 km s"^ for H/3, 380±38 km 
s-i for Hei A6678, 390±23 km s'^ for Hei A7065, 480±27 km s'^ for Hen A4686 and 490±66 km s"^ for Niii AA 4641,4642. 
The velocity separations of Ha and H/3 are the same within the uncertainties of the measurements, and so are those of Hex 
A 6678 and Hei A 7065. The similarity of the velocity separations for Hen A 4686 and Niii A A 4641,4642 implies that these lines 
are emitted at a similar distance from the compact object, consistent with the theory of Bowen fluorescence, which requires 
Hen for the pumping mechanism (see McClintock, Canizares & Tarter 1975; Kallman & McCray 1980; Schachter, Filippenko 
& Kahn 1989). Thus, we have observed a clear trend in which the velocity separations of the double-peaked lines increase 
with line ionisation, from the neutral hydrogen lines (fs 250 km s~^), to the neutral helium lines (w 400 km s~^) and finally 
to the higher-ionisation (i.e. ionised helium and Bowen Niii) lines (~ 500 km s^^). 

The optical continuum of GX 339—4 was bright during our 1998 observations, suggesting the presence of an active, 
bright accretion disk. If the double-peaked lines are emitted from the accretion disk, as commonly interpreted, we can use the 
peak separations to infer the radial distance of the emission regions from the compact object. Moreover, as higher-ionisation 
lines are emitted from hotter gases, we can also use the peak separations to map the surface temperature distribution of the 
accretion disk. 

Assuming that the lines are emitted from a geometrically thin, Keplerian accretion disk, we can use the line velocities 
to constrain the system parameters. In Figure 2 we plot the expected velocity separation as a function of distance from the 
compact object for systems with a compact object of mass 4.0, 7.0 and 10.0 M© and for various orbital inclinations. We also 
mark in the figure the expected orbital separation of the systems for various orbital periods, assuming that the companion 
star is a Roche-lobe-filling main-sequence star. As the size of the accretion disk cannot exceed the Roche lobe of the compact 
star, which is smaller than the orbital separation, we can easily see that the orbital inclination of GX 339—4 is well below 30° 
if the orbital period is 14.8 hr (Callanan et al. 1992). Moreover, the mass of the compact object is unlikely to exceed 10 M©. 
However, the orbital inclination and the mass of the compact object may be higher if the orbital period is longer. If we adopt 
the BHC average mass of 7 M,:? (see e.g. McClintock 1998) as the mass of the compact object in GX 339—4 and assume that 
the Balmer lines are emitted near the rim of the accretion disk and that the orbital period is 14.8 hr, then we deduce an 
orbital inclination of about 15° for the system. 

We note that no large velocity variations of the centres of the disk lines were detected in our 1998 April and August 
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Figure 2. The theoretical velocity sepciration of the line peaks as a function of the distance from a compact star of mass equal to 4.0, 7.0 
and 10.0 Mq (from left to right), assuming a geometrically thin, Keplerian accretion disk. The curves correspond, from top to bottom, 
to the orbital inclinations i = 75, 60, 45, 30 and 15° respectively. The horizontal lines are the velocity separations deduced from our 1998 
observations. The short vertical lines mark the orbital separation of the system for multiples of the suspected 14.8-hr orbital period. 



observations. This may be further evidence of a low orbital inclination. Additional evidence is provided by the flat radio 
spectrum (Fender etal. 1997), which we interpret as optically thick emission from jets oriented close to the line of sight. 



4 SPECTRUM IN THE LOW-HARD STATE 
4.1 General features 

Figure 3 shows the high-resolution (1.3 A) spectra that we obtained on 1999 April 12. For comparison we also show the 

lower-resolution (3 A) summed spectrum from the 1997 May observations, when the system was in a previous low-hard state. 
In the 1999 spectra both the Hq and the H/3 lines arc prominent, with equivalent widths (EWs) of —8.5 ±0.3 and —4.0 ±0.5 A 
respectively. The Ho lino appears to be single-peaked, in contrast to the more clearly double-peaked profile observed in the 
high-soft state. For H/J, the signal-to-noise of the spectrum does not allow us to determine the line profile. However, we do 
not see strong evidence of two peaks. 

The other prominent line is Hen A 4686 (EW — —6.0 ± 0.5 A). It shows two clearly resolved peaks, unlike the Balmor 
lines. The Niii AA 4641,4642 lines are detected (with a limit to the EW > —1.8 A) but are much weaker than in 1998 April 
and August (EW ~ —(2.5 — 3.0) A) during the high-soft state. Other lines that can be idcntifiod arc weak Hei A 4922, and 
He I A 6678 emission. In the 1997 May spectrum, the Hei A 6678 line is also found, and the Nii A 6505 line is visible. 



4.2 Are the emission lines single-peaked or double-peaked during the low-hard state? 

Medium-resolution (« 2.5 A) spectra obtained by Smith, Filippenko & Leonard (1999) during the low-hard state in 1996 
May showed a broad, flat-topped Ha line. The profile resembles a line with two unresolved peaks, and Smith, Filippenko & 
Leonard (1999) interpreted the Ha line as being double-peaked. However, our 1997 May observations, when GX 339—4 was 
still in its low-hard state, did not show a double-peaked Ha line. Because of the low resolution (3 A) of our spectra, we were 
unable to verify whether the Ha line was single- or doubled-peaked. The higher-resolution (1.3 A) spectrum that we obtained 
on 1999 April 12 allows us to carry out a more quantitative analysis. 

The profiles of the Ha and the He ii A 4686 lines in the high-soft and low-hard states are shown in Figure 4. The Ha 
profiles of GRO J1655— 40 obtained during its high-soft state in 1996 June and the hard X-ray outburst in 1994 September 
are also shown for comparison. 

We have used the normalised difference spectra to examine the differences in the flux distribution in the line profiles from 
each observing epoch. For each spectrum, we removed the continuum by fitting and subtracting a low-order polynomial in 
the velocity range between —4000 km s~^ and ±4000 km s~^. The residual line flux was then normalised to unity. Finally, 
the line spectrum obtained in the low-haxd state was subtracted from the spectrum obtained in the 1998 high-soft state. 

In the left panel of Figure 5, we show the difference profile of the Ha lines for the 1998/1999 observations. The two peaks 
in the 1998 observations arc clearly seen above the noise level. The peaks are also visible, though less prominently, in the Ha 
difference profile for the 1997/1998 observations (Fig. 5, middle panel). This suggests that the Ha profiles are different in the 
high-soft and the low-hard states. While the high-soft state is characterised by a double-peaked profile, the line probably has 
a single-peaked profile in the low-hard state. 
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The Hen A 4686 line is doublc-pcaked in both our 1998 and 1999 spectra, making the difference profile (Fig. 5, right 
panel) less straightforward to interpret. The overall shape results from the difference between the widths of the lines. 

Our observations show that both the Ha and the He ii A 4686 lines are double-peaked in the high-soft state. In the low-hard 
state, the Hen A 4686 line clearly shows two peaks, but the profile of the Ha line is better described as single-peaked. This 
indicates (i) different emission regions for Ha in the high-soft and the low-hard states, (ii) different line-formation mechanisms 
in the two states, and (iii) different response of the outer accretion disk under irradiation by soft and hard X-rays. These 
issues will be discussed in more detail in the later sections. 

4.3 Comments on the line profiles and line widths 

4.3.1 Hen \4686 and Nui X\ 4641,4642 

In spite of the similar double-peaked profiles of He 11 A 4686 in the two X-ray spectral states, the peak separations are clearly 
diflferent: 490±66 km s"'^ for the high-soft state and 600±50 km s~^ for the low-hard state (Fig. 4, top-right panel). The line 
widths are also different in the two states (Fig. 5, right panel). The larger separation implies that the emission region is at a 
smaller radial distance from the compact object. The larger peak separation in the low-hard state may be due to the fact that 
the soft X-ray component was weak and unable to set up the temperature-inversion layer (see §5) in the outer disk regions. 

We note that the Niii AA 4641,4642 fines in 1999 April (Fig. 3) were much weaker than in 1998 (see Fig. 3 in Paper I). 
The Niii AA 4641,4642 lines are due to Bowen fluorescence, which is pumped by the On: A 374 multiplet resulting from the 
cascade O III 2p3d ('^Pl^') 2p3p (^S,^P,^D) 2p3s (''^Po) 2p^ (-^Pz) transitions primarily pumped by the UV Hen A 303.78 
(Lyman-a) line (see e.g. Deguchi 1985; Schachter, Filippenko & Kahn 1989). The Hen Lyman-a line can be produced by 
the photoionisation-recombination process, and is therefore associated with the He 11 A 4686 line. Hence, one might expect to 
observe strong Niii AA 4641,4642 emission when Hen A 4686 is prominent. This is what we have observed in GX 339—4 and 
other BHCs (e.g. GRO J1655-40 and RXTE J1550-564) during the high-soft states. The situation is different during the 
low-hard state of GX 339—4 in 1999 April: the Niii AA 4641,4642 lines are weak, yet the Hen A 4686 remains prominent. A 
similar phenomenon was also observed in the BHC GRO J0422-I-32 during a hard X-ray outburst (Casares et al. 1995). 

It was suggested that the absence of Niii A A 4641,4642 lines in GRO J0422-I-32 was due to insufficient X-ray irradiation 
(Casares et al. 1995). One would therefore expect both the He 11 A 4686 and the Niii AA 4641,4642 lines to be weak. Here (and 
in various other systems), our observations show a strong Hen A 4686 line both in the high-soft and low-hard states, but 
significantly weaker Niii AA 4641,4642 lines in the low-hard state. Clearly, there must also be other factors. We propose that 
the weakening of the Niii AA 4641,4642 emission from BHCs in the low-hard state is due to optical depth effects. When the 
Hen region becomes optically thick to UV continuum, the Hen Lyman-a and Oiii A 374 photons arc quenched. If it is still 
relatively transparent to the optical continuum. Hen A 4686 photons can escape. Without efficient pumping by Hen Lyman-a 
and On: A 374 photons, the Nni AA 4641,4642 lines are suppressed. If this is the situation in GX 339—4, the photoionisation- 
recombination processes which produce the He 11 A 4686 lines may be different during the high-soft and low- hard X-ray states 
(cf. Case A, Case B and Case C recombination processes, see Osterbrock 1989; Xu et al. 1992). 

4-3.2 Line width of the Ha emission line 

The profiles of the Ha and the He 11 A4686 lines observed from GX 339—4 in the high-soft and low-hard states are shown in 
the first three panels of Figure 4. The width of Ha is similar in the two states, in contrast to Hen A4686, in which case the 
line width increases with the hardness of the X-rays. 

We interpret the double-peaked lines in the high-soft state of GX 339—4 (and in the high-soft state of other BHCs, 
like GRO J1655— 40) as originating from the temperature-inversion layer created by soft X-ray irradiation on the accretion 
disk surface. The single-peaked (flat-topped or round-topped) Balmer lines observed in the low-hard states of GX 339—4 are 
probably emitted from an outflow, possibly a dense wind from the evaporating atmosphere of the accretion disk. Theoretical 
studies (Murray et al. 1995; Murray & Chiang 1998) have shown that if the outflow/wind has a substantial velocity component 
parallel to the disk plane, single-peaked lines with widths similar to the Keplerian disk velocities can be produced. The profiles 
will appear to be round-topped or flat-topped, depending on the outflow velocity and density profiles (see Chapter 14 in Mihalas 
1978), as in the case of emission lines from windy massive stars. 

Single-peaked Balmer emission lines were also observed from GRO J1655— 40 in 1994 August - September, in between 
strong hard X-ray flares. (Soria, Wu & Hunstead 2000). On that occasion, the lines were much narrower than the double- 
peaked lines observed during the high-soft state of the system (Fig. 4, bottom-right panel). The narrow, single-peaked lines from 
GRO J1655— 40 are probably emission from an extended, optically thin cocoon above the disk plane, with lower rotational 
velocity than the disk or the disk wind. In the high-inclination system GRO J 1655— 40 (i ~ 70°), most of this extended, 
optically-thin envelope was seen projected onto the sky and was therefore a source of emission lines. Narrow emission lines 
have not been observed in the low-inclination system GX 339—4, where the accretion disk is seen almost face-on and any thin 
cocoon would be seen projected onto the disk surface. 



5 LINE EMISSION REGION 
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5.1 A plane-parallel model 

We consider a simple plane-parallel model (Fig. 6) to illustrate that a temperature-inversion layer can be set up at the surface 
of an accretion disk under strong soft X-ray illumination but not under hard X-ray illumination. The model is a modification 
of the model constructed by Milne (1926) for the illumination of a star by its companion in a close binary system. In Milne's 
original formulation the incident radiation is at optical wavelengths, and the cooling is via emission of optical radiation. In 
BHCs, the incident radiation on the accretion disk consists of soft and hard X-rays emitted from regions close to the central 
compact object. We therefore need to take this difference into account. 

In our model, we consider that the incident radiation has soft and hard X-ray components, for which the opacities are 
respectively larger and smaller than for the optical radiation. This modification is essential in the study of accretion disks 
irradiatively heated by X-rays, as soft X-rays are easily absorbed by neutral and weakly ionised matter (via bound-free 
transitions) at the disk surface, while hard X-rays are attenuated only at large depths where the matter density is sufficiently 
high. The other assumptions are the same as those in the Milne (1926) model, namely (i) the emission region has a plane- 
parallel geometry, (ii) the system is in radiative equilibrium, and (iii) the principle of superposition is applicable. 

As the radiative-transfer equations are linear, the local temperature stratification is given by 



T(r) = 



Bx(r)-|-Bd(T) 



1/4 



-,1/4 



-Mr) 



(1) 



(see Milne 1926), where B^{t) is the component of the radiation due to irradiative heating, Sd(T) is the component due to 

viscous heating in the disk in the absence of X-ray irradiation, r is the optical depth, and o is the Stofan-Boltzmann constant. 
We assume that the disk is cooled by emission of optical radiation, and therefore, t refers to the optical depth in the optical 
band, along the direction normal to the disk surface. The disk component is given by 



3 



2Tto 



(2) 



(e.g. Dubus et al. 1999), where the viscous energy flux of the disk Sd = aT^g/Tv and rtot is the total opacity of the accretion 
disk in the vertical direction. The effective temperature Tg^, which is determined by balancing the energy generated by viscous 
heating and by radiative loss (see e.g. Smak 1984), is given roughly by 



efl 



(3) 



where v is the viscosity coefficient, E is the surface density of the accretion disk, and Qk is the Keplerian angular frequency. 

Suppose the incident radiation consists of parallel beams of soft and hard X-rays, with effective fluxes nSs and irSh 
per unit area normal to the beams, and making an angle a with the normal to the disk plane. The beams are absorbed 
exponentially. We denote the absorption coefficient of the soft X-rays by fe/t, and that of the hard X-rays by khK, where k is 
the absorption coefficient of the optical radiation. 

At the optical depth t, the soft and hard X-ray fluxes are attenuated to ttSbC^''"^^^''"' cosa and yrShe"*'"'^'''''^" cosa 
respectively. The absorbed flux of X-rays in depth dr is therefore -!T{kBSBe~''^'^'""'" + khShe~''^'^ dr. Radiative equihbrium 
requires 



+ feh^he 



')+[ 

J 47r 



dn I(t, h) = 47rBx(r) 



(4) 



where I{t,h) is the intensity of the optical radiation from the irradiatively heated disk at the depth t in the direction 
(= cos~^ /i). 

Recall the radiative-transfer equation 



t^-^I{r,tx) = 7(T,M)-i?x(T) . (5) 

Multiplying the above equation by a differential solid angle dfl, integrating with respect to r and using the condition of 

radiative equilibrium, we have 



[ dnnI{T,n)= n cos a (S, e''"'^ " -h Sh e-*!-" ") 



(6) 



With appropriate boundary conditions, the radiative-transfer and radiative-equilibrium equations can be solved by the 
method of successive approximations (Milne 1921, 1926). In the limit of a semi-infinite plane the first-approximation solution 
is 



Bx(r) 



cos a ^ 1 
kg 2 



-kgT sec cx 
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, _ , , cos a 1 



COS Oi 1 \ —/gj^T sec OL 



(7) 



The derivation of the above equation will be presented elsewhere (Wu & Soria, in preparation). 

In the second approximation, Bx(t) is assumed to have the same functional form as that of the first-approximation 
solution, i.e. 



B^{r) = a-K 6"*="="'=" - 6h 6"'=''"="'^" , 



(8) 



where a, 6b and 6h are constants to be determined by the boundary conditions. For a semi-infinite slab opaque at optical 
wavelengths, the emergent radiation is the Laplace transform of Bx(r), i.e. 



7(0, m) 



= lim 
= a — 6s 





cos a 



-6h 



kh 



kh 



The condition of no inward component of I{t, /i) at the boundary (t = 0) implies 

dp 1(0, n) = 



/: 

and 

J dfj, fj, 7(0, /i) = . 

With these boundary conditions, we obtain from the radiative-equilibrium condition 



(9) 



(10) 



(11) 



a - 6s 



, cos a 1 , , ^ 
2 - I — — 1 ln(l + fes sec a) 



, cos a 1 , /, , N 
2- ( — — ) In(l + A:hseca) 



Similarly, the radiative-transfer equation is integrated, yielding 



o — 2 6s 



ln(l + ks sec a) 



-2 6h 



- (fesS's + fehSh) 



ln(l + kh sec a) 



= (^Sa + Sh) COS a 



(12) 



(13) 



Using the trivial conditions 6s — » when Ss — > and 6h — > when Sh — > 0, we can solve the radiative-equilibrium and 
radiative-transfer equations (12) and (13) for the constants a, 6s and 6h, and obtain 



6s — ^ 



6h = -khSh 



kh 



Ma) , 



and 



fe5s l^]Ma) + khSh i ^ ) fh{a) 



where the two functions /s(a) and /h(tt) are given by 



(14) 
(15) 

(16) 



Ma) = 
and 

Ma) = 



, cos a \ / cos a \ / cos a 1 , , 



, cosa \ / cosa \ / cos a 1 i , /, , x 



(17) 



(18) 



We now define a hardness parameter ^ = Sh/Ss and a total X-ray flux S:ii_ = Sg + Sh- Then, we have 



2 



-ksT sec a 



+A;h/h(a) 



1 + ^ 



-fchT sec a 



(19) 
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5.2 The temperature-inversion layer 

Using the plane-parallel model, we first investigate the temperature profile, as a function of optical depth r, in the absence 
of radiation due to viscous dissipation in the accretion disk. In this case the temperature profile is determined purely by the 
soft and hard X-ray irradiation. Cold disk matter has a larger opacity for soft X-rays than for optical radiation because of 
bound-free absorption. Hard X-rays are not easily absorbed, and electron scattering is the main source of hard X-ray opacity. 
Wo assume a simplified situation in which the an average opacity is assigned for each waveband, and consider the parameters 
ks — 5.0 and fch = 0.01 for the soft and hard components respectively. 

We have found that a strong temperature-inversion layer can be set up at the accretion-disk surface when the incident 
angle of the X-rays approaches grazing incidence (i.e. a — + 90°). For normal incidence, the temperature inversion is weak and 
unimportant (Fig. 7a) . The temperature inversion is caused mainly by soft X-ray irradiation, but the more penetrating hard 
X-rays tend to heat up deeper layers (Fig. 7b). 

Next we investigate the situation where viscous dissipation is also present, i.e. the temperature profile is determined 
jointly by viscous heating in the accretion disk {Sd > 0) and by X-ray irradiation (Sx > 0). We consider a particular situation 
in which the total flux resulting from the two processes is constant but the relative strength of the two processes varies. We 
show in Figure 7c the temperature profile for irradiation from a very soft X-ray source (with a hardness parameter ^ = 0.01) 
at an incident angle a — 87°. When the X-ray component dominates (i.e. F = S^/Sd > 1), a temperature- inversion layer is set 
up. When the X-rays are weak the temperature profile is determined by viscous heating. Figure 7d shows the corresponding 
temperature profile for the case in which the X-rays are hard (denoted by a hardness parameter ^ = 100). The temperature 
profiles are similar in spite of the different relative strength of the viscous heating and irradiation heating. When the incident 
angle is smaller (say a = 60°, see Fig. 7e), the effects due to irradiation become even less important. Although hard X-ray 
irradiation cannot significantly alter the overall temperature profile, it does increase the temperature uniformly at all optical 
depths — a direct consequence of the small opacity for hard X-rays (feh -C 1) and of our approximation that feh is depth 
independent. It is worth noting that the apparent effects at very large optical depths (t » 1, see Fig. 7e) when the hard X-ray 
fiux is much stronger than the viscous dissipation fiux arise from the assumptions of a semi-infinite plane and single-sided 
illumination. The situation would be different for a finite slab with both the top and the bottom surfaces being illuminated; 
this model will be discussed elsewhere (Wu & Soria, in preparation). 

We finally consider a situation in which the total X-ray flux is flxed at 10 times the viscous dissipation flux but its 
hardness is allowed to vary. This corresponds to strong X-ray irradiation with a constant bolometric X-ray luminosity and 
with spectral evolution. Our model shows that the temperature- inversion layer becomes more prominent as the X-ray hardness 
parameter ^ decreases (Fig. 7f). Also, while the temperature at small optical depths (near the accretion-disk surface) has 
changed substantially, the temperature at large optical depths (in the mid-disk plane) is only modified slightly. 

In summary, the vertical temperature profile of an accretion disk is sensitive to soft X-ray irradiation if the X-ray flux is 
stronger than the viscous dissipation flux, but it is not significantly altered by hard X-ray irradiation. Soft X-ray irradiation 
results in the formation of a temperature-inversion layer on the accretion-disk surface; hard X-rays tend to provide uniform 
heating at all depths in the disk. 



5.3 Line formation 

The transfer equation for non-polarised radiation in a semi-infinite plane-parallel atmosphere, with complete redistribution 
and consideration of non-LTE, is 

H-J-I,^ = —{k^ -\- K-i^<pv)Iv -\- K^S^ -\- Ki^(puS^ (20) 

dz 

(see Mihalas 1978). The subscripts "c" and "l" denote the continuum and the line-centre respectively, k,^ and are therefore 
the opacities of the continuum and line respectively. Similarly, and are the continuum and line-centre source functions. 
(j)v is the line profile function, and Iv is the intensity at frequency v. 

Without loss of generality, we simply let /i = 1 (0 = 0°). As an approximation we can rewrite the transfer equation at 

the line-centre frequency as 

= -(Kc +ic +iL , (21) 

where the emissivities and j-^ are related to the corresponding source functions as usual. In the neighbourhood of the line, 
4>v ^ Q, and therefore we have the transfer equation for the continuum 

^Jc = -Kc-^c+ic (22) 
(e.g. Tucker 1976). It follows 

Jo Jo 

where lo is the background intensity, and the total optical depth r is the sum of the optical depths of the line and the 
continuum, i.e. t = + ''l • Generally, ^ . Hence, t w . A line appears in emission if > Ici and in absorption if 
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7l < Iq.I^ a, homogeneous medium, this is equivalent to a criterion > lo + Iq for hue emission and < Io + Iq for hne 
absorption, where 



For systems not far away from LTE, the criterion above can be approximated by B{T) > lo + 1^ for hne emission, and 
B{T) < lo + Ic for hne absorption. (Here, T is the thermal temperature of the region where the line originates, and B{T) is 
the local Planck function.) An inhomogeneous medium can be considered as a stack of homogeneous strata, and the criterion 
above is still applicable. 

As a qualitative illustration, we may consider a simple model in which the accretion disk consists of two strata: a central 

disk plane and a surface layer. Now, the source function of the lino from the surface layer is denoted by S-^ , and the intensity 
of the continuum emission by I^. The Planck function of the surface layer is B{T). The emission from the disk beneath 
can be considered as the background 7o. For emission lines to be observed, the line- formation region is either an extended 
atmosphere/corona transparent to the continuum or a temperature-inversion layer above the surface of a geometrically thin 
and opaque accretion disk. 

A corona can be considered as a geometrically extended surface layer. For certain viewing geometries (e.g. when the 
disk is viewed edge-on), the background disk emission Jo is unimportant. As the corona is transparent to the continuum 
emission, its intensity 7^ is smaller than the local Planck function B{T). The line emission is, however, opaque, implying that 
= B{T) > lo -\- If,. Hence, the line-emission criterion is satisfied. 

If the surface layer is geometrically thin and the disk is opaque, then lo cannot be neglected. Since the temperature- 
inversion layer has higher temperatures than the disk, B{T) > Jo. Moreover, as the temperature-inversion layer is optically 
thin to the continuum, B{T) 3> I^ - For LTE, = B{T), and the criterion B{T) > lo + 1^ is therefore satisfied, provided 
that the temperature of the surface layer is significantly higher than the temperature of the disk. For non-LTE, the line 
centre has a brightness temperature much higher than the local thermal temperature T of the surface layer, and hence 



Wc now consider the temperature profiles of a region on a disk under moderate and strong irradiation of soft X-rays 
with various incident angles, and make use of them to investigate whether or not emission lines will be formed. In the case 
of moderate irradiation, we set F to be 0.42 (Fig. 8, left panel). Temperature inversion is insignificant in this case, especially 
when the X-ray incident angle a ^ 45°, and the temperatures at the disk mid-plane are higher than the temperatures of 
the surface layers. For LTE the intensity of the lines from the surface layer is roughly the local Planck function of the layer. 
As the intensity of the emission from the disk lo is roughly given by the effective temperature of the disk, there will not be 
emission lines in this case. If the line emission process is non-LTE, then the brightness temperature of the line is higher than 
the thermal temperatures of the surface layer and the disk, and hence it is possible to form emission lines. 

In the case of strong X-ray irradiation (F — 1000) the surface layer is significantly hotter than the disk (Fig. 8, right 
panel). Temperature inversion is stronger for larger incident angles. As the Planck function of the temperature-inversion layer 
is greater than the Planck function of the disk where the optical depth is equal to unity, the lines are in emission. 



5.4 Line-emission region of GX 339—4 

In our simple plane-parallel model we have assumed parametric X-ray opacities that scale with the optical opacity, and 
considered a particular case, in which feg = 5 and kh = 0.01. We now examine whether or not these values are appropriate for 
GX 339—4. The effective X-ray opacity is given by Kx.eff = [/«x(«;x -|- Kes)]^''^ (Rybicki & Lightman 1979). The X-ray opacity 
is Kx « 47.9 (0.8 keV/Ef cm^g-^ for 0.8 keV < E < 8.0 keV and Kx « 0.17 (8.0 keV/Ef cm^g~^ for E > 8.0 keV (e.g. 
Sincell & Krolik 1997). The opacity due to electron scattering is Kes = 0.4 cm^g~^. At T ^ 3 x 10^-10'* K the optical opacity 
is mainly due to atomic bound-bound transitions. An exact treatment requires the details of atomic processes and also the 
chemical composition of the accreting matter (see e.g. CoUin-Souffrin & Dumont 1990; Alexander, Johnson & Rypma 1983). 
As an approximation, we consider the Rosseland mean opacity k « 0.7 (p/lO^^g cm^^)(r/10^K)~^'®cm^g~^ given in Prank, 
King & Raine (1992), and assume that it can be used to describe the opacity of the optical/UV radiation from the accretion 
disk. The soft component of the X-rays from GX 339—4 is a thermal black body with E ~ 1-3 keV. If we take i5 ~ 1.5 keV, 
we obtain a soft X-ray opacity fta ~ 7 cm^g~^. The hard component is a power law extending between i5 ~ 10 100 kcV. If we 
take a value E ^ 40 keV, we have a hard X-ray opacity Kh ~ 0.02 cm^g~^. The optical opacity is ~ 1 cm^g~^. Our choices 
of fcs and fch are therefore appropriate. 

In a previous section, we have proposed that the double-peaked lines seen in GX 339—4 during the high-soft state originate 
from a temperature-inversion layer above an irradiatively heated, Keplerian-like accretion disk. Now, we are going to use the 
plane-parallel model to test whether this scenario is consistent with our 1998 observations. 

For black-hole binaries, most of the X-rays are emitted from a small region within 1000 Schwarzschild radii from the 
central black hole. The linear size of the X-ray emission region is therefore < 10® cm. For a system with a period greater than 
12 hr, the radius of the accretion disk would extend beyond ~ lO^^cm, and hence the incident angle of the X-rays on the 
accretion disk is close to 90° . 

In our calculation, we consider a standard ShakurarSunyaev disk around a 7-M0 black hole under irradiation of X-rays 
at an incident angle of 87°. We assume a viscosity parameter of avis = 1 for the accretion disk, and the X-ray luminosity 
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from the central (point) source given by Lx = rjMc^ , where M is the mass-accretion rate and c the speed of hght. We adopt 
an efficiency parameter 77 = 0.1 following Prank, King & Raino (1992), and a mass-accretion rate of 5 x 10^^ g s~^. Taking 
into account that most of the soft X-rays are reflected at grazing incidence, a reflective efficiency of 99% (see de Jong, van 
Paradijs & Augustcijn 1996) is assumed. 

In the left panel of Figure 9 we show the temperature structure of the accretion disk illuminated by soft X-rays with a 
hardness ratio ^ = 0.02. Strong temperature inversion occurs at radii R ^ 10^ cm, with surface temperatures (at r = 0) about 
10'' - 10^ K at Ji ~ 10" cm. The effective temperature (evaluated at t w 2/3) is about a factor of 2 — 3 lower. We also show 
the locations of the line-emission regions deduced from the peak separations of the lines during the 1998 high-soft state, for 
an assumed orbital inclination of i = 15° (cf. Fig. 2). The temperature distribution resulted from our model calculation is in 
agreement with the temperatures of the formation regions of doublc-pcakcd hydrogen and helium lines. 

The middle panel of Figure 9 shows the same accretion disk but under hard X-ray irradiation. The hardness ratio is 
^ = 50. The soft component of the X-rays is weak, and is insufficient to create a prominent temperature inversion. For LTE 
processes, the lines omitted from the disk should be suppressed. It is worth noting that although the composition of the 
X-rays is very different in the two cases, the effective temperature of the disk is similar. While the emission lines are formed in 
the temperature-inversion layer, the continuum emission is characterised by the effective temperature. One therefore expects 
that the optical brightness in the two cases will not differ significantly. Thus, our model shows that the optical continuum is 
relatively insensitive to the X-ray hardness, as observed in GX 339—4, provided that other factors do not alter significantly 
when the X-ray hardness varies. 

For comparison we also consider an accretion disk irradiated by moderately hard {(, = 10) X-rays and assume that only 
50% of the X-rays are reflected (Fig. 9, right panel). As shown, the entire accretion disk is severely heated, because the energy 
flux deposited by the X-rays is larger than the energy flux generated by viscous heating in the disk. There is a temperature 
inversion at the disk surface because of a non-negligible soft component, but it is not as strong as the case where the X- 
rays have an overwhelmingly strong soft component. Under such severe irradiativo heating by the hard X-rays, the structure 
of the accretion disk should deviate significantly from the structures deduced from the Shakura-Sunyaev prescription. The 
radiative-transfer and hydrodynamics equations must be solved simultaneously in order to obtain a self-consistent temperature 
structure of the accretion disk. (This is beyond the scope of the present paper.) Nevertheless, we have shown here that the 
accretion disk responds very differently to X-ray irradiation of different hardness. 

We suspect that a prolonged low-hard state causes parts of the accretion disk in GX 339—4 to be heated up. The outer disk 
inflates to become an extended isothermal atmosphere/corona, and the temperature- inversion layer is quenched in this region. 
The X-ray heating may also drive a disk wind. The Balmer lines, which axe emitted from the extended atmosphere/corona or 
wind, lose their Keplerian-disk signature and become single-peaked (cf. Murray etal. 1995). The Hen A 4686 line is emitted 
at a smaller radial distance from the compact object, where the accreting material is still relatively confined to the disk plane, 
and the residual soft component of the X-rays causes a local temperature inversion. Hen A 4686 may therefore retain its 
double-peaked profile. 

Although the plane-parallel model considered here produces results consistent with our spectroscopic observations of 

GX 339—4 in the 1998 high-soft state, we are aware of its limitations. In particular, non-LTE effects are not considered. 
Moreover, we have assumed a "semi-grey" approach, and the atomic processes that determine the radiative equilibrium and 
govern the radiative transfer are not treated explicitly. The plane-parallel model is sufficient for the purpose of a qualitative 
illustration, but a proper study should treat all these effects in more detail. 

Double-peaked lines are also observed from some black-hole candidates in quiescence, for example A0620-00. However, 
we suspect that the mechanism of omission is very different in the two states. We showed in §5.3 that line emission can occur 
in two physical conditions: either from a temperature-inversion layer on the surface of optically thick gas, or from optically 
thin gas; in the former case, the lines are superimposed on the contimiuni emission. In the case of GX339-4 in outburst, the 
double-peaked optical emission lines come from the irradiated surface of an optically thick disk. In the case of A0620— 00 
and similar systems in quiescence (or more generally at very low accretion rate and X-ray luminosity), the accretion disk is 
probably optically thin, hence it is hot. In this case, omission lines do not require a temperature-inversion layer to be formed, 
hence they are not caused by soft X-ray irradiation. It is worth noting that GX 339—4 in quiescence (Orosz 2000, private 
communication) also shows Balmer emission while the V brightness is at least five magnitudes lower than in the soft and hard 
X-ray states, suggesting that the accretion disk is optically thin in the continuum. 



6 CONCLUSIONS 

Our observations have shown that the profiles of the emission lines from GX 339—4 are different in the high-soft and low-hard 
X-ray spectral states. In particular, the hydrogen Balmer lines have single-peaked profiles in the low-hard state but double- 
peaked profiles in the high-soft state. The Hen A 4686 line is doublod-peaked in both states. Our interpretation is that soft 
and hard X-rays have different heating effects on the accretion disk and, as a result, the line-formation mechanisms in the 
two X-ray states are different. We have constructed a simple plane-parallel model to illustrate the differences in the two cases. 
Our model has shown that a strong temperature-inversion layer can be formed when the X-rays arc soft and they illuminate 
the disk at an angle approaching grazing incidence. If the double-peaked lines that we observed during the high-soft state 
are emitted from the temperature-inversion layer on a geometrically thin, opaque, Keplerian accretion disk, we deduce that 
GX 339—4 is a low-inclination system. The orbital inclination will be about 15° , for an orbital period of 14.8 hours. 
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Figure 3. The high-resolution (1.3 A) spectra that we obtained on 1999 April 12 when the system was in a low-hard state. About five 
days later the system entered an X-ray quiescent state. The spectra are centred at Ha (top panel) and at H/3/Heii (middle panel). The 
lower-resolution (3 A) summed spectra, centred at Ha, obtained in 1997 May (low-hard state), is also shown for comparison (bottom 
panel). 
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Figure 4. Profiles of Ha (top-left panel) and Hen A 4686 (top-right panel) during the high-soft state in 1998 April and the low-hard 
state in 1999 April (represented by solid and dotted lines respectively). We smoothed the 1999 low- hard state profiles by averaging 
over three data points. We show the comparison of the Ha-line profiles for the 1997 and 1999 low-hard states (dotted and solid lines 
respectively, bottom-left panel). The 1997 data are not smoothed. Profiles of the Ha line from the high-orbital-inclination (i 70°) BHC 

GRO .11655—40 arc shown (bottom-right panel). The broad, double-peaked Ha profile of GRO J1655— 40 (solid line) was obtained on 
1996 June 10 when the system was in a high-soft state; the narrow Ha profile (dotted line) was obtained on 1994 September 3 when it 
was active in the hard X-ray bands. The normalisations of the continua in the GRO J1655— 40 spectra fall below unity because of the 
broad Ha absorption component (Soria, Wu & Hunstead 2000). 
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Figure 5. Normalised difference Ha line profiles of the 1998/1999 and 1997/1998 observations (left and middle panels respectively). The 
locations of the peaks of the line obtained in the 1998 observations are marked (vertical lines) . We calculated the diflference profiles by 
subtracting the profile obtained in the 1997/1999 low-hard states (Si) from those obtained in the 1998 high-soft state (S2). The integrated 
line fluxes (after continuum removal) are normalised to unity. The 1-<t noise level (dotted lines) is calculated from the continuum in the 
velocity ranges between —4000 km s~^ and —1000 km s"'^, and between 1000 km s~^and 4000 km s~^. The difference profile of He 11 
A 4686 for the 1998/1999 observations is also shown (right panel). Note that the excess in the flux at velocities < —1500 km is due 
to the Nlll AA 4641,4642 lines, which were more prominent in 1998 than in 1999. 




Figure 6. The geometry of the modifled IVIilne (1921) plane-parallel model. 
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Figure 7. (a) Temperature profile (T cx B^/'^) as a function of optical depth t for X-ray incident angles a = 87° (solid line), 60° (dotted 
line), 30° (dash-dotted line) and 0° (dashed line). The viscous energy flux from the disk = 0, the X-ray incident flux Sx = 1, the 
X-ray hardness parameter ^ = 1, and the total optical depth rtot = 20. The two scaling parameters for the X-ray opacities feg = 5.0 
and fch = 0.01. (b) Temperature as a function of optical depth t for X-ray hardness parameters ^ = 0.01 (solid line), 0.1 (dotted line), 

1 (dash-dotted line), 10 (dashed line) and 100 (dash-dot-dot-dottcd line). The incident angle a is fixed at 87°. The other parameters 
are the same as in (a), (c) Temperature as a function of optical depth r for non-zero viscous disk flux S^. The curves correspond to 
illumination-strength parameters F = S^/S^ = 100 (solid line), 10 (dotted line), 1 (dash-dotted line), 0.1 (dashed line) and 0.01 (dash- 
dot-dot-dottcd line) with t; = 0.01 (i.e. strong soft X-ray illumination). We fix 5x -I- = 1.0 and a = 87°. The other parameters are the 
same as in the previous cases, (d) Same as (c) but with !; = 100 (i.e. strong hard X-ray illumination), (e) Same as (d) but with a lower 
incident angle a = 60°. (f) Same as (b) except that here we have F = 10 (i.e. > 0). 
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Figure 8. Temperature profile as a function of optical depth r for moderate (F = Sx/Sj = 0.42) and strong (F = 1000) soft X-ray 
irradiation (left and right panels respectively). The X-ray hardness parameter is ^ = 0.01, the two scaling parameters for the X-ray 
opacities are kg = 5.0 and feh = 0.01. 5x + Sd is fixed to l.The X-ray incident angles are a = 0°, 45° and 87° (lines from top to bottom). 
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Figure 9. The temperature structures of an irradiativcly heated accretion disk calculated from our plane-parallel model. The surface 
temperature (at t = 0) of the disk is represented by the dotted line; the effective temperature (at t = 2/3) by the thick solid line; 
and the temperature at t = 10 by the dashed line. We have assumed a Shakura-Sunyaev disk with a viscosity parameter Ovis = 1- 
For comparison we also show the effective temperature (thin solid line) of the same Shalcura-Sunyaev disk without X-ray heating. The 
mass of the black hole is 7 Mq; the mass-accretion rate is 5 x 10^^ g s~^; and the incident angle of the X-rays is fixed to be 87°. In 
the left panel, the disk is irradiated by soft X-rays with a hardness ratio ^ = 0.02, and 99% of the X-rays are reflected. We also mark 
the radii at which the Balmcr, He I A 6678, He I A 7065, He II A 4686 and Nlll AA 4641,4642 lines arc emitted, inferred from the peak 
separation observed from GX 339—4 in 1998. In determining the location of the linc-cmission region, we assume an orbital inclination 
i of 15°, a constraint imposed by the 14.8-hr period (§3 and Callanan ctal. 1992), and that the lines arc emitted from a geometrically 
thin, Keplerian accretion disk (see Fig. 2). If the period is longer or the orbital inclination is larger, the linc-cmission regions will be at 
some larger radii (Sec Fig. 2). Balmer lines are emitted at temperatures ~ (1 — 2) X 10'* K; He I lines at ~ (2 — 4) X lO** K; and He II lines 
at ~ (5 - 10) X 10-* K (see e.g. Cox & Tucker 1969; Dopita & Sutherland 1996; Osterbrock 1989). The Bowen Niii AA 4641,4642 lines 
are formed in regions where He ii A 304 photons are abundant. As shown, the temperatures at the inversion layer are consistent with the 
temperatures at which these lines are emitted. In the middle panel, the accretion disk is irradiated by X-rays with a hardness ratio of 
^ = 50, and 99% of the X-rays are reflected. In the right panel, the hardness ratio ^ = 10, and only 50% of the X-rays are reflected. 
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